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Abstract 

The role of the Weibel instability is investigated for the first time in the context of the large-scale 
magnetic reconnection problem. A late-time evolution of magnetic reconnection in relativistic pair 
plasmas is demonstrated by particle-in-cell (PIC) simulations. In the outflow regions, powerful 
reconnection jet piles up the magnetic fields and then a tangential discontinuity appears there. 
Further downstream, it is found that the two-dimensional extension of the relativistic Weibel 
instability generates electro-magnetic fields, which are comparable to the anti-parallel or piled- 
up fields. In a microscopic viewpoint, the instability allows plasma's multiple interactions with 
the discontinuity. In a macroscopic viewpoint, the instability leads to rapid expansion of the 
current sheet and then the reconnection jet front further propagates into the downstream. Possible 
application to the three-dimensional case is briefly discussed. 
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I. INTRODUCTION 

Magnetic reconnection is widely recognized as a fundamental physical mechanism in 
collisionless plasmas. Consuming the magnetic field energy in the inflow region, it releases 
the energy to the kinetic energy of plasma particles. It is an effective engine for magnetic 
dissipation, plasma heating, or particle acceleration. Recently, the relativistic extension of 
magnetic reconnection has received attention for its role in various high-energy astrophysical 
places — active galactic nuclei [H [2j [3], pulsars jU El [6], gamma ray bursts [7J |8], and 
magnetars [HI EI] • The mechanism of relativistic reconnection still remains unclear as well as 
the conventional non-relativistic reconnection, but recent kinetic simulations start to reveal 
its features. It is demonstrated that relativistic pair plasma reconnection is a powerful 
particle accelerator [TH [TJl [131 HH US]; an d that it keeps fast reconnection rate [161 E] 
despite the lack of the Hall physics [18]. In principle, magnetic reconnection is a relatively 
large scale process — the typical speed of reconnection jet is the Alfven velocity Va or the 
light speed ~ c, and the typical time scale is several tens of characteristic time scale; X/Va 
or A/c, where A is the typical spatial scale of the field reversal. 

On the other hand, in the context of gamma ray bursts [T9l 120] or extra galactic jets, 
the Weibel-type two-stream instability in relativistic plasmas has attracted recent attention, 
too. The Weibel instability [21] is an electromagnetic instability, that arises from plasma 
anisotropy. Since it quickly generates magnetic fields, it is a likely origin of magnetic field in 
the synchrotron source near collisionless shocks or near relativistic jet fronts. Series of PIC 
simulations [221 E31 EH ES, EEl [2TJ, EH] successfully demonstrate the magnetic generation via 
jet penetration or plasma shell collision in a weakly or nonmagnetized plasma. The Weibel 
magnetic structure evolves into long-durated filament-like magnetic structure, whose energy 
is approximately 10% of equi-partition energy. In-situ particle acceleration is also reported, 
but its detailed mechanism still remains unclear [221 ESI EH]- On the theoretical side, the 
conventional Weibel instability, which propagates into the transverse direction from plasma 
anisotropy, has been extended to relativistic temperatures [3TJ [321 EH]- Meanwhile, its two- 
dimensional extension, the electromagnetic counter streaming instability, has been studied 
in relativistic counter streaming conditions [221 EH EQ] . In general, the Weibel instabilities 
are microscale phenomena, whose scales are characterized by the plasma frequency u p = 
[Aline 2 /m] 1 / 2 ; uj~ 1 in time and by (ck/u p ) in space. 
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In the context of magnetic reconnection, since magnetic reconnection expels powerful 
outflow jets from the reconnecting region, it is quite possible that the jets interact with 
pre-existing plasmas, and then excite an anisotropy-driven instability. In fact, Daughton 
& Karimabadi [34| reported generation of out-of-plane magnetic field via firehose type in- 
stability in their non-relativistic pair plasma reconnection, although its role in reconnection 
remains unclear. 

In the present paper, we study the role of the Weibel instability in the reconnection con- 
text. We carry out two-dimensional PIC simulations of relativistic pair plasma reconnection, 
and we find that the relativistic counter-streaming Weibel instability generates out-of-plane 
magnetic fields in the downstream region of reconnection outflow. We discuss the properties 
of the instability, and then we investigate how the Weibel instability affects micro- and macro 
physics of magnetic reconnection. The paper consists of the following sections. In section II 
we describe our simulation setup. In section III we present the two-dimensional simulation 
results, and then in section IV we investigate the properties of the instability. In section V 
we discuss how the Weibel instability effects plasma dynamics, both in microscopic particle 
motion and macroscopic reconnection structure. The last section contains discussion and 
the summary. 

II. SIMULATION SETUP 

We carry out two-dimensional particle-in-cell (PIC) simulations in a current sheet con- 
figuration. As an initial condition, we employ a relativistic extension of the Harris model 
in GSM-like geometry. The magnetic field, plasma density and plasma distribution func- 
tions are described by B = B tanh(z/A)x, d(z) = d cosh~ 2 (z/A) = (jpn ) cosh~ 2 (z/ A) and 
f s oc d(z) exp[—jp{E — P s u y } /T]. In the above equations, B is the magnitude of antiparallel 
magnetic field, A is the typical thickness of the current sheet, do is the lab- frame number den- 
sity in the current sheet, n is the proper number density, the subscript s denotes the species 
('p' for positrons, 'e' for electrons), (3 P = —j3 e = j3 is the dimensionless drift velocity, 7/3 is 
the Lorentz factor for (3 (7^ = [1 — /3 2 ] -1 / 2 ), e is the particle energy, u is the relativistic four 
velocity of u = [1 — (V/c) 2 ]™ 1 / 2 • v and T is the proper temperature including the Boltzmann 
constant. We set T = mc 2 and j3 = 0.3, respectively. In addition, a uniform background 
plasma is added to the system in order to supply plasmas in the reconnection inflow region. 
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Its number density and temperature are db g /do = 5% and Tb g /mc 2 = 0.1, respectively. In 
general, the velocity of the reconnection outflow jet is known to be approximately the Alfven 
velocity in the inflow region. In the present study, we choose a low-density (5%) background 
population to obtain fast outflow. Notice that the Harris model with uniform background 
plasmas exactly satisfies an equilibrium. 

The system consists of 1568 (a;) x 768 (z) grids and the typical scale of the current sheet A 
is set to 10 grids. Since we consider periodic boundaries in the x direction, and since there 
are two current layers in the periodic z direction, the boundaries of the main simulation 
domain are located at x — ±76. 8A and z = ±19. 2A. We use 7.5 x 10 7 super particles in this 
simulation. One cell contains 6.6 x 10 2 particles at the center of the current sheet. During 
the very early stage we impose a small artificial electric field E y around (x, z) = (0, ±3A). 
The resultant E x B flow compresses the current sheet, and then it triggers reconnection 
around the center of the main simulation domain. The typical spatial ranges of the trigger 
field E y are set to (Ax,Az) ~ (±2A,±A). The trigger field soon vanishes after t/r c = 
(10 — 15), where r c = A/c is the light transit time. We discuss the physics of reconnection in 
the late stage of t/r c = (60 — 120), which is not influenced by this initial perturbation. These 
conditions are similar to the author's previous study [H] ; but we use a larger simulation box 
to discuss late-time structure without boundary effects. We call this reference run 'run A'. 
The total energy is conserved within an error of 0.1% throughout the simulation run, after 
the initial trigger force vanishes. 

III. RESULTS 

Due to the initial perturbation, magnetic reconnection takes place around the center of 
the main simulation domain. Magnetic field lines start to reconnect at t/r c ~ 50, and outflow 
jets into the ±x directions appear. Figure [l]a show a snapshot of the right half of the main 
simulation domain at t/r c = 80. The left half (—76.8 < x/X < 0) is not presented. The 
reconnection outflow is well developed at this stage, and its speed is up to ~ 0.7c. Along 
the neutral line, magnetic fields are piled-up in front of the dense plasma region of the 
current sheet around x/X ~ 18. Its peak amplitude is B z /B ~ 1.5, and there is a relatively 
sharp boundary between the pileup magnetic field and the pre-existing dense plasma in 
the downstream. This boundary is a tangential discontinuity (hereafter TD in short), and 
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we discuss "upstream" and "downstream" based on the TD throughout this paper. The 
propagation speed of the TD (Vtd ~ 0.65c) is slightly slower than the average velocity 
of local plasmas. The typical plasma density in the simulation frame is d/do ~ 0.1 — 0.2 
in the upstream, d/do ~ 2.5 at the downstream side of the TD and then it decreases to 
d/do ~ 1 in the further downstream region. Figure shows the out-of-plane magnetic field 
(By) structure and the relevant current system in the reconnection outflow region, which 
is indicated by the rectangle in Figure [T]a. The characteristic structure of B y is observed, 
and its maximum amplitude is B y ~ 0.6B . In Figure [T]c we observe charge separation at 
the same place and the vertical E z structure (Hereafter the term "vertical" means the z 
direction). The E z explains both the motional field for B y and the electrostatic field by 
the charge separation. The time development of the B y structure along the neutral plane 
is presented in Figure [2j These B y fields suddenly appear after t/r c ~ 64 and then they 
exponentially grow until they saturate after t/r c ~ 80. The instability looks like a convective 
mode, traveling into the +x-direction. However, actually, it is nearly non-convective purely- 
growing mode in the frame of the plasma average flow. The linear growth rate measured 
by B y growth is T c Ui ~ 1.7-1.8 x 10 _1 or Ui/Q p ~ 5.2-5.5 x 10~ 2 , where Q p is the typical 
plasma frequency in the system. The typical spatial scales are 7A-10A (x) and ~ 2 A (z). 
Careful observation show that the instability has a two-dimensional rectangular structure. 
In Figure [T]6, we see the weak negative regions on the upper side, on the lower side and on the 
right side [colored in light blue; (x, z) ~ (24, ±1.5), (30, 0) in unit of A] of the characteristic 
positive region [inner orange region; ~ (24,0)]. Similarly, weak positive regions are located 
in the vicinity of the characteristic negative region [inner blue region; ~ (21,0)]. We find 
that these structures are generated by the two-dimensional Weibel-type instability. In this 
case, plasmas are highly anisotropic along the x-direction, mainly because the TD pushes 
away the pre-existing plasmas, and because the reconnection outflow jet penetrates into this 
region. Therefore, the situation is similar to jet injection [2U or relativistic counter- 
stream [22], and magnetic generation near the shock [20] in pair plasmas. The instability 
resides inside the current sheet, where the plasma frequency is high. In addition, the Weibel 
instability prefers an unmagnetized region, and so an inner current sheet is an ideal place 
for the instability. The current structure and the charge separation structure indicates the 
nature of the Weibel-type activity. As schematically explained in Medvedev & Loeb [20], 
small By fluctuation leads to the ^-displacement of ±x-streaming plasmas, and then the 
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resultant x-current structure ±5 J x continues to enhance 8B y . Thus, the Weibel instability 
generates magnetic field which is perpendicular to the direction of the anisotropy, and then 
it leads to the reduction of anisotropy. 

Figure [3] shows the plasma distribution function of 1.4 x 10 5 particles in the Weibel active 
region (22 < x/X < 26, -2 < z/X < 2) at t/r c = 80. The left part of the distribution 
function is almost identical to the initial distribution of pre-existing plasmas, but the right 
part are highly elongated due to both reflected plasmas and upstream-origin accelerated 
particles. The average plasma velocity is (v x ) ~ 0.49c, where ( ) means the average value 
among the particles. On the contrary, the plasma fluid velocity (a Lorentz transformation 
velocity to the rest frame, where the plasma momentum flow is zero) is ~ 0.65c ~ Vtd- The 
two velocities differ due to highly asymmetric plasma distribution. An integrated plasma 
temperature is as follows; (T x ,T y ,T' z ) = ((mu' x v' x ), (mu' y v' y ), (mu' z v' z )) = mc 2 (2.3, 1.0, 1.5), 
where T", mu' and v* are the temperature, momentum and velocity in the rest frame of 
plasmas. After the reconnection jets start from the X-type region, the plasma anisotropy in 
the downstream region grows in time, until the Weibel instability appears. The anisotropy 
stays at the same level after the instability occurs, because x-momentum is continually 
supplied from the upstream side. 

Snapshots of the field properties along the neutral line; the pileup field B z {bold line), 
the reconnection electric field E y {dashed line), and the Weibel magnetic field B y [thin line) 
are presented in Figures [6]a-<i Note that the Weibel fields are observed in the local frame of 
plasma average flow. They propagate to the +x direction and its speed is slightly slower than 
the speed of the TD [Vtd)- Therefore, sometimes the Weibel fields are caught up by the TD. 
For example, the positive B y region around x/X ~ 24 at t/r c = 80 is nearly caught by the 
TD around x/X ~ 34 at t/r c = 100. At the same time, new Weibel fields are continuously 
generated in the further downstream region. We discuss the late time development of the 
Weibel fields later. The subpartition of Weibel field energy (B y + E x + E z )/8n to the local 
plasma kinetic energy saturates around 8-12% in run A. 

IV. LINEAR ANALYSIS 

In order to study the properties of the instability, we have solved the dispersion relation 
by linearizing relativistic fluid equations. In the Weibel region, plasmas mainly consist of 
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three different components; (i) pre-existing plasmas in the Harris current sheet, (ii) current 
sheet plasmas, which are reflected by the TD, and (iii) upstream-origin plasmas, which 
are originally from the reconnection inflow region. Since we set a low plasma density in the 
reconnection inflow region (5% of the Harris current sheet), the third population is relatively 
smaller than the other two. Therefore, we employ the counter-streaming model of (i) and (ii) 
in order to evaluate the Weibel instability. We extend Kazimura et al. [22] 's fluid theory for 
counter-streaming four fluids (streaming/counter-streaming positrons and electrons), which 
was originally developed by Califano et al. [29J. Although Kazimura et al. |22j ignored 
the plasma pressure effect, it is here included. We employ the following relativistic fluid 
equations; 

"^f (Psa + e sa ) ( 7^ +V sa ■ V) V^ a = -Vpsa + lsaqsaTisa ^+^x5j-^ (lsaq S an sa E ■ V sa + ) , 

(1) 

where p is isotropic plasma pressure, e is the fluid internal energy, the subscript a denotes 
two kind of streams ('1' for streaming fluids, and '2' for counter-streaming fluids), and 
7 sa = [1 — (v sa / 'c) 2 ] _1//2 is the relevant Lorentz factor. We also use the continuity equation 
and Maxwell equations 

d 

Qlilsansa) + V • (j sa n aa v aa ) = 0, (2) 
„ - 4vr v-^ 1 dE 

V X B = y y ^saqsaUsaVsa H 7TT (3) 

c ^— ' c at 

s=e,p a=l,2 

We assume the adiabatic gas condition in order to close the equation 

p sa oc n r sa , e sa = n sa mc 2 + - - - p sa (5) 

where T = 5/3-4/3 is the polytropic index of adiabatic gas, We consider a two dimensional 
perturbation 5f oc 5f exp(ik x x+ik z z —iut) , where k = (k x , k z ) is the wavevector and u is the 
complex frequency, and then we linearize all equations for four fluids. Then, we numerically 
calculate the growth rate (Im u) for arbitrary wavevector k = (k x , k z ) by solving a matrix 
problem. For simplicity, the following assumptions are used; 

v p i = v el = V T d, v p2 = v e2 = -V T d (6) 
n p i = n p2 = n el = n e2 = n (7) 

Ppl = P P 2 = Pel = Pe2 = UqIJIC 2 . (8) 



We assume that plasma density is homogeneous, two counter-streams are symmetric, the 
frame is set to the co- moving frame of the TD, considering that the TD completely re- 
flects the momentum of pre-existing plasmas. In the present case, the simulation data 
shows k = (k x ,k z ) ~ (ti; p /c) (0.15, 0.75) in the frame of interest. Because of the complex- 
ity in the simulation system, this analysis does not exactly describe the instability. The 
density gradient of plasmas, the current sheet thickness, the wavelength of the instability 
are all comparable, the Weibel region moves to the x-direction slightly slower than the TD 
(Vtd ~ 0.65c), the local average velocity and the local fluid velocity differs, and local plasma 
velocities depend on the distance from the TD. However, the goal of our simple theory is to 
roughly understand the physics. 

Figure [4] shows the dispersion relation of the two-dimensional mode for k z = 5k x . The 
linear analysis (bold line) and the simulation data are in good agreement. Further investiga- 
tion shows that the maximum growth rate is on an order of 0.01-0.02u;p with the relativistic 
temperature of T = mc 2 , and that the cut-off (decline of the growth rate) is rather sensitive 
to the counter-streaming velocity. The obtained mode is purely growing, and it has an elec- 
tromagnetic feature. Because of the mathematical symmetry, we obtain the other oblique 
modes for (±k x , ±k z ) with the same growth rates. Therefore, the two-dimensional rectan- 
gular structure is obtained by superimposing these oblique modes. The change separation 
structure (Figure [T]c) in the simulation frame can be explained by the ^-displacement by 
the instability. It reflects both the density gradient inside the current sheet and the Lorentz 
boost of the fast outflow streams. The electrostatic component of the instability is relatively 
small. 

We can also obtain the growth rate of the instability in counter-streaming cold beams 
[22] by dropping the plasma pressure effect. (One can remove pressure-related terms from 
eq. [I] and employ e sa = n sa mc 2 instead of eqs. [HJ) For comparison, the growth rate of the 
cold-beam limit is also presented in Figure [4] (dashed line). Obviously, the instability grows 
substantially slower than the cold-beam limit. One interpretation is that imposing plasma 
pressure means the reduction of the anisotropy. In a high temperature limit where the four 
velocity of the counter-streams is relatively negligible, the distribution becomes close to a 
single isotropic distribution. Another interpretation is that the relativistic pressure effect 
slows down the growth rate, as discussed in the relativistic studies on the one-dimensional 
Weibel instability [3TJ E2J ES]- In relativistic temperature regime, it is known that the Weibel 
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instability is re-scaled by Im uj < uj p /^ 1 / 2 in time and {^ l / 2 ck)/uj p in space, where 7 is the 
typical Lorentz factor of plasma maximum energy [HI]. From the viewpoint of relativistic 
fluids, the enthalpy term in equation [T] increases an effective inertia, and then it slows 
down the growth rate of the instability. In the present case, the term yields {^ sa / c) 2 {p sa + 
e sa ) ~ ( nm )lsa(^ + [r/(r — 1)] [T sa / (mc 2 )}) ~ 9n sa m. Since it replaces the mass term inside 
the plasma frequency, the instability in a relativistic hot plasma grows slower than the 
instability in the cold beam limit by a factor of \/9 ~ 3. By comparing the enthalpy term in 
relativistically hot limit (~ 4p/c 2 ) and in cold-beam limit (~ nm), we obtain the slow-down 
factor of p 1 / 2 . This is consistent with the scaling of the one-dimensional Weibel instability, 
by a factor of 7 1 / 2 . In summary, the counter-streaming Weibel-type instability slows down 
by the inertia effect of relativistic pressure. Roughly speaking, the instability is similarly 
re-scaled by a factor of 7 1 ' 2 , as the one-dimensional Weibel instability. 

The panels in Figure [5] present growth rates of the obtained unstable modes as a function 
of k = (k x , k z ). Both relativistic pressure case (Fig. [HJa) and the cold-beam limit (Fig. ^p) 
are shown. The one-dimensional mode along k x = is the conventional Weibel instability, 
which has electromagnetic features. On the other hand, the mode along k z = is the 
electrostatic counter-streaming instability. The typical mode in our simulation is (k x , k z ) = 
(0.15,0.75) with some amount of ambiguity. As seen in Figure [5]a), the obtained mode 
is rather close to the one-dimensional Weibel instability. It is important to note that the 
oblique mode grows slightly faster than the one-dimensional Weibel instability, and this is 
a signature of the counter-streaming Weibel-type instability. The central region and the 
right half of Figure [5]a are mainly occupied by the electrostatic-like mode. Their growth 
rate is even faster, however, since our theory depends on the isotropic fluid pressure and 
the adiabatic condition (eq. [5]), we think that our theory may be invalid, especially in the 
short wavelength range of (c/u> p )\k\ > 1. In addition, in the high-pressure regime, the two 
counter-streaming distributions overlap each other. All these conditions are unfavorable to 
describe the electrostatic modes in the parallel direction. Meanwhile, the cold beam limit 
(Fig. ^p) seems to be in good agreement with the Saito & Sakai [SOj's work, based on 
Kazimura et al. [22] 's theory. Their counter-streaming velocity 0.5c is comparable to ours 
of 0.65c. 



9 



V. EFFECT OF THE WEIBEL INSTABILITY 



In this section, we investigate how the Weibel instability affects the micro dynamics 
of plasmas motion and global dynamics of reconnection. First we focus on the plasma 
motion near/in the downstream region, because the Weibel instability occurs only in the 
downstream side of the TD. Two characteristic regions will affect plasma motion — the TD 
and the Weibel region. Before the Weibel instability appears, particles in the downstream 
region are meandering in the current sheet. Once they are hit or reflected by the TD, 
they constantly travel into the +x-direction because ^-momentum is conserved in the Harris 
current sheet configuration without B y . 

So, what happens after the Weibel magnetic fields B y appears? In order to study plasma 
motion around the two characteristic regions, we select 10 5 super particles (5 x 10 4 pairs) 
that satisfies the following conditions; they are (i) found in the piled- up region (6 < x/X < 
10,-2 < z/X < 2) at t/r c = 60 and (ii) found in the Weibel region (22 < x/X < 26,-2 < 
z/X < 2) at t/r c = 80. The x-ranges of these regions are indicated by arrows in Figures 
|6]a, b. Then, the spatial distribution of the selected particles are investigated. The panels in 
Figure [7] show the distribution of the selected positrons at t/r c = 100. The right panel shows 
the distribution of fast positrons. The relevant rc-range is indicated by the dashed arrow in 
Figure |6]c. The left panel presents slow positrons, whose x- velocity v x is slower than that of 
the TD; vtd — 0.65c. The panel contains 10 4 slow positrons (20% of positrons). 

Roughly speaking, these positrons can be classified into the following three groups. The 
first group is moving-away positrons, who travels faster than the TD to the +x direction. The 
rectangle (indicated by the white dashed line in Figure [7]&) is an approximate location of the 
selected particles at t/r c = 100. After they are hit by the TD, they continue to escape into 
the +x region, faster than Vtd- The second group is found around 33 < x/X < 40, z/X ~ ±2 
in both two panels in Figure [7| They are located along the magnetic field line, which are 
connected to the TD. They have relatively small population, and they do not always escape 
into the +x direction. We discuss the field-line modulation and the current sheet expansion 
later in this section. The last group is found along the neutral line (z ~ 0) in Figure [7]a. 
These particles are affected by the Weibel instability. Their ^-locations are positive around 
40 < x/X < 44, and negative around 36 < x/X < 40. These ^-displacements are due to the 
Weibel instability; the effect of the out-of-plane field B y (Refer to Fig. lolc for the polarity 
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of By). The high density region near the TD (x/X < 36, z/X ~ 2 in Figure [7]a) contains 
both the second class of positrons along the field lines and the third class of Weibel-affected 
positrons. They are soon reflected by the TD, and then we see their reflection in the other 
high density region near the TD (x/X < 3Q,z/X ~ —2 in Figure [7]&). We note that the 
magnetic field near the TD is not vertical, but rather tilted into +y direction, because the 
TD hit the positive B y region at this time. Since 20% of selected positrons are in Figure [7^, 
the third group has relatively large population. 

We pick up 1.7 x 10 3 positrons from them, which are found in the vicinity of the TD 
(x/X < 33, —2 < z/X < 2) at t/r c = 100, and then we examine their trajectories. The 
selection mainly consists of the third class of Weibel- modulated positrons, because the second 
class has fewer population. Properties of two typical trajectories are shown in Figure [6]e ; / 
as a function of x. We call them positron A (solid line) and positron B (dotted line). 
Throughout the simulation period (0 < t/r c < 120), they stay in the narrow region of 
—2 < z/X < 2. In Figure [6]e, marks show the particle position at the selected time stages 
for comparison with Figures ^p,-d. Particle A starts from x/X ~ 33 in the — x direction. Its 
energy is originally e ~ 3mc 2 . Around t/r c ~ 60, it collides with the TD and then it turns 
its way to the +x direction. Its energy increases to e ~ 5mc 2 via the interaction with the 
TD. Figure [6J7 show the field properties B y (solid line) and E y (dashed line) at positron 
A's position. Near x/X ~ 22, it feels relatively strong B y in the Weibel-active region, and 
then its x-momentum is transported to z momentum (also, to y— momentum through the 
meandering motion). Consequently, its x- velocity starts to slow down, as indicated by the 
arrow in Figure [6]f. Since the TD travels relatively fast (vtd ~ 0.65c), the TD eventually 
catches up positron A and hits it again. This time, the positron A gains more energy by the 
motional electric field E y , because more magnetic fields are piled up than during the first 
impact. Now its energy goes up to e/mc 2 ~ 8. So, the Weibel field slow down the escaping 
particles, and then it enables multiple interactions with the TD. We even find three-time or 
four-time interactions with the TD in selected positron trajectories. For example, positron 
B (dotted line) are hit by the TD twice (x/X ~ 31,44) after the first interactions with the 
TD at x/X ~ 9. If we study further long-time evolution, these positrons will be hit by the 
TD multiple times. 

Interestingly, it seems that only low-energy particles are reflected by the Weibel region; 
high-energy particles are insensitive to the Weibel fields and then they easily escape to the 
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+x direction, when their energy e/mc 2 exceeds 7-8. The typical kinetic energy gain by the 
TD reflection ranges up to e/mc 2 ~ 10-20, due to enhanced pileup electric field. Therefore, 
among the reflected particles, low-energy particles are trapped between the TD and the 
Weibel region, and then they are heated by the multiple interactions with the TD. They can 
not escape into the outflow region, until they become energetic enough — their gyro radii 
exceeds the scale of the Weibel structure {^c/uj c ^> rf^c/ujp). So, this result indicates that 
the Weibel instability enhances plasma heating in the downstream side of the TD, while it is 
not likely to enhance high-energy particle acceleration. The threshold energy will increase as 
the system condition becomes more relativistic, because the electron skin depth is relatively 
larger by default, and because the typical scale of the Weibel instability becomes even larger 
by a factor of 7 1 / 2 or p 1 ^ 2 . 

Next, in order to study the Weibel mode effect to global reconnection structure, we carried 
out another simulation. The new run (run B) starts from intermediate data of run A at 
t/r c = 60, and then we artificially reduce B y to in run B. The quadrupolar magnetic fields 
near the reconnecting region [221 EZ] do not appear in pair plasma reconnection, because 
Hall physics depends on different ion and electron masses. In antiparallel configurations, 
magnetic reconnection involves B x and B z , and only the Weibel instability or anisotropy- 
driven instabilities generates B y . Because of the artificial reduction of B y , the system slightly 
loses energy. The total energy in run B is smaller than the total energy in run A by 0.15% 
at the end of the simulation (t/r c = 120). This energy loss is much smaller than typical 
kinetic, magnetic, and internal energies. The panels in Figure [8] show the late-time evolution 
of the two runs at t/r c = 120. For comparison, the upper halves present the results of run A, 
while the lower halves present those of run B. At this time, there should be no substantial 
boundary effects, since the evolution time of the reconnection region is smaller than a wave 
propagation time to the ^-periodic boundaries. In addition, the magnetic field lines at the 
A-type region come from z/X ~ ±15 around x/X = ±76.8 as seen in Figure |8]a. The 
reconnection still involves magnetic flux inside the main simulation domain. 

One can see the difference in global structure in Figure [8]a. The current sheet seems 
to be broadened in run A, while it remains thin in run B. The position of the TD front 
differs, too. The TD is located at x/X ~ 46 in run A. On the contrary, the TD is located 
at x/ A ~ 43 in run B — the TD front can not penetrate into the x-direction as run A. In 
Figures [6]a-d, the B z profile along the neutral line in run B is also presented by a dotted 
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line. The peak plasma density at the downstream of the TD is d/do ~ 2.4 in run A, while 
plasmas are much compressed near the TD; d/d Q ~ 5 in run B. These differences can be 
explained by the current sheet expansion by the Weibel instability. In run A, the Weibel 
instability transfers some of plasma x-momentum into ^-momentum (and also ^/-momentum 
via meandering motion). Therefore, it reduces plasma x-pressure in the downstream region, 
and the TD front can move further distance into x-direction. The increased ^-momentum 
leads to the current sheet expansion, and then magnetic field lines become more round. The 
current sheet continues to expand as long as plasma ^-momentum is continuously supplied 
through the Weibel activity, from reconnection outflow in the upstream region. On the 
contrary, plasmas and anti-parallel field lines are rather confined near the current sheet 
in run B. The TD can not penetrate into the outflow direction as run A, therefore, the 
maximum amplitude of the pileup field is stronger (See Figs. ^jp,-d for the B z profiles in run 
B) and the inflow speed near the TD is slightly slower than in run A (Fig. $d). 




The total reconnected flux £|-B 2 | along the neutral plane is the same in runs A and B, 
because it is controlled by the physics of the upstream region; magnetic reconnection near 
the X points. Around (x, z) ~ (0, 0), we observe a small magnetic loop in the current sheet. 
This is not a projection of the X-point, but a secondary magnetic island, which appears after 
t/r c ~ 105. Its formation mechanism is unclear. Since these secondary islands are also found 
in non-relativistic reconnection in pair plasmas [34] , the island formation will be common 
feature in a low-density current sheet in large scale simulations. Anyway, we do not see 
noticeable difference — both run A and run B are almost same in the upstream region. At 
this stage the Weibel activity is not likely affecting the central reconnection region, because 
information can not have widely spread over the central reconnection region. For example, 
an information of large -B y -perturbation at x/X ~ 20 at t/r c = 80 (Fig. [2]) can not arrive 
at the X-type region before t/r c = 100. On a longer time scale, the Weibel instability may 
have an impact on the reconnection rate since it changes the magnetic field line topology by 
expanding the current sheet |40j . 

In the case of run A, Figure [8^ and Figure |8]c show the downstream field structure in 
more detail. In Figure[8]fr, the magnetic field lines near x/\ ~ ±1 at the right boundaries are 
set to connect to x/X = ±1 at periodic boundaries. The field line shifts to z/X ~ 3 at the 
thickest point due to the current sheet expansion in run A, while the field line stays around 
z/X ~ — 1 in run B. In this stage, the Weibel instability is also active outside the neutral line 
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around zj A ~ 2, as well as along the neutral line. The xz current system is well developed 
around the B y regions. Charge distribution (Fig. [8]c) is correlated to the a;-current system 
(Fig. |8]6); we see positron-rich J x > region and electron- rich J x < region. Compared 
with the early stage in Figure [TJ these structures are rather elongated into x-direction. This 
is consistent with many studies on Weibel instability; elongated "filament" structure or 
current channels are commonly observed in well-developed stage of the Weibel instability. 
In run B, there is no current system in the xz plane. Regarding the outflow structure in 
Figure [8]ci, plasma flow is rather bifurcated in run A, due to the ^-displacement of plasmas. 
We can see a significant difference in the J y current structure in Figure |8]e. In run A, the 
current region is located in front of the broadened plasma region. The energy conversion 
J ■ E mainly takes place in the vertical current front there. In run B, the ^/-current structure 
is enhanced around the small spot near the TD, and then energy conversion takes place 
there. 

Regarding the composition of the energy in the system of interest, two runs slightly dif- 
fer in accordance with the field line topology; the summary of the upstream reconnection 
field energy £(-B 2 + B 2 + Ey)/8ir is almost same in both two runs, but run A has more 
field energy (120%) in £(-B 2 + E 2 )/8n compared with run B. The total plasma kinetic en- 
ergy £(7 — l)mc 2 is almost same, however, in run A, a slight percentage of them (0.5%) 
are converted again into the Weibel-related field energy; £(-B 2 + E 2 + E 2 )/8n. The total 
amount of Weibel-related field energy is equivalent to ~ 5.3A 2 (£> 2 /87r). This is substantially 
smaller (10 _1 -10 -2 ) than that of the upstream-related field energies, because the Weibel 
active region is relatively small. 

Figure [9] presents energy spectra in the regions of interest. The spectra of two runs look 
similar, too. However, in order to distinguish the difference clearer, these spectra are divided 
into two parts by the TD; in the downstream region of the TD, and in the upstream region 
of the TD. Note that the TD is located in the further downstream in run A. We observe a 
high-energy nonthermal tails in their spectra in the upstream side. This is due to dc particle 
acceleration or piled-up acceleration in the upstream side [HJ [H] . In the mid-energy range 
(20 < e/mc 2 < 60), run A has slightly more high-energy population. We think this is due 
to the larger volume of the upstream region. Since particles can stay longer inside the larger 
upstream region, or the main site of particle acceleration, more particles are accelerated 
into high energy range in run A. In the low-energy range around e/mc 2 ~ 10, We expected 
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that plasma heating is enhanced through multiple interaction by the TD in run A, but the 
enhancement is too small (even in linear scaling). One reason is that the Weibel region is 
too small. Furthermore, the Weibel region not only reflects the escaping particles, but also 
it hits the pre-existing particles. So, as a result, the net effect will be small. Meanwhile, in 
run B, plasma population is slightly enhanced around e/mc 2 ~ 20. It is difficult to discuss 
this energy range, because too many effects are relevant. 

VI. DISCUSSIONS AND SUMMARY 

In ion-electron reconnection, it is well known that quadropolar out-of-plane fields B y 
appears in the vicinity of the A-type region [351 ESI EI]- However, quadropolar structure 
disappears in pair plasmas with an equal temperature [TO]. In the present case, the out- 
of-plane fields are found in the downstream of the reconnection outflow region. Daughton 
& Karimabadi [34| reported similar structure in non-relativistic pair plasma reconnection, 
and they argued that it comes from some type of firehose instability. Since both the fire- 
hose instability and the Weibel instability belong to the anisotropy-driven instabilities, the 
generation of the out-of-plane fields in the downstream would be common feature in pair 
plasma reconnection. We identified that the out-of-plane field B y is generated by the Weibel 
instability downstream of the TD. A comparison with another run (with B y artificially sup- 
pressed) demonstrates that the Weibel instability leads to a significant modulation of the 
downstream structure; further penetration of outflows, the current sheet expansion and the 
bifurcated downstream jets. We expect that the formation of "T-shaped current sheet" [12] 
can be explained by the current sheet expansion by the Weibel instability. 

In three dimensions, the Weibel instability also generates the vertical magnetic field ±5B Z 
and the out-of-plane electric field ±5E y . The Weibel instability will lead to a filament-like 
development of the x-currents, involving small-scale reconnection of perturbed magnetic field 
lines [38] ■ The plasma's x-momentum will be transferred to y-momentum as well as to z- 
momentum. We expect that the TD penetrates further downstream into the outflow region, 
because plasma x-pressure will be more efficiently scattered. Meanwhile, the current sheet 
expansion in the downstream region may be less apparent, because all of the x-momentum 
is not transformed into ^-momentum. Regarding the particle acceleration, we will observe 
more high-energy particles, because the TD will further penetrate into downstream and then 
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the upstream acceleration site expands. In addition, high-energy particles from the upstream 
region may also be affected by the Weibel fields. In two dimensional case, such high energy 
particles are rather insensitive to the Weibel region, and low-energy reflected particles bounce 
between the TD and the Weibel-active region. However, in three dimensions, the Weibel 
magnetic field can affect high energy particles, especially when it is antiparallel (B z < 0) to 
the pile-up field. Thus, some higher-energy particles may bounce between the TD and the 
Weibel-active region. Along with the expansion of the upstream acceleration site, particle 
acceleration is likely to be enhanced. 

Furthermore, we should consider all other instabilities in three dimensions. It is known 
that the relativistic drift kink instability (RDKI) quickly modulates the current sheet in a 
relativistic pair plasmas [TH 139] . Its typical growth rate in this configuration is t c uj.i ~ 0.1, 
while T c Ui ~ 0.03 for the tearing instability. Although the RDKI grows slower than the 
Weibel instability, the RDKI is a macro instability, and it may inhibit the reconnection 
process by modulating the current sheet [T3] , while the Weibel instability is the sub-product 
of the reconnection outflow. Since the RDKI slowly widens the current sheet, unmagnetized 
or weakly magnetized region becomes wider. Therefore, we expect that the Weibel instability 
is active in a wider region inside the modulated current sheet. In addition, since the Weibel 
instability involves y-structure, repeated collision between the TD and the Weibel fields may 
lead to the instability of the TD front in the xy plane (e.g. the interchange instability of the 
reconnection jet front [H]). The Weibel instability in three dimensions will be an interesting 
problem to challenge. 

The Weibel instability will also occur under non-relativistic temperature condition of 
T <C mc 2 . In this regime, usually the electron skin depth becomes smaller than the other 
scales like electron gyro radius, and so the Weibel instability occurs in a shorter time/spatial 
scale in reconnection. As long as it occurs inside reconnection outflow structure, the physics 
will be the same. 

In ion-electron plasmas, the Weibel instability will work for electrons, and then it may 
contribute to quick electron heating. Although it is not clear whether the sharp TD is 
formed in the outflow region in ion-electron plasmas, the multiple interaction with the TD 
will also be possible. Similarly, enhanced heating may also occur near the fast shock or the 
other discontinuities. In solar cases, it is reported that hard X-ray emission comes from the 
small spot near the loop top of magnetic field lines (12] . Shock- related electron heating may 
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take place in the downstream of reconnection outflow, and the Weibel instability and the 
relevant bounce effect may contribute to the enhanced electron heating there. In addition, 
the Weibel instability may play a role in quick electron heating inside the nanoflare jets, in 
the context of coronal heating problem. 

In more generalized configuration of magnetic reconnection, for example, in magnetic re- 
connection with uniform guide field B y , the situation will differ substantially. Since the guide 
field By scatters x-momentum into z-momentum, the wavevector of the Weibel instability 
is likely to be in the y direction. However, since the outflow is slower than the antiparallel 
case, and since the ambient magnetic field B y stabilizes the instability, the Weibel instability 
will be less active. The situation will be more complicated in relativistic pair plasmas, be- 
cause charge neutrality often breaks down in the outflow region [15]. Therefore, how plasma 
anisotropy disappears in the guide field case remains to be solved. 

On the viewpoint of energetics, the ultimate energy source of the Weibel instability is 
the plasma bulk energy of the reconnection jet, which is expelled by the magnetic energy in 
the inflow region. Initial magnetic energy is converted to plasma energy of reconnection jet, 
and partially to magnetic energy of the Weibel region. Then, the Weibel activity modifies 
the downstream reconnection structure, which potentially changes the downstream energy 
conversion process once again! 

Finally, let us briefly summarize this paper. We investigated the role of the Weibel 
instability in the reconnection context. We demonstrated the following new results; (1) the 
Weibel instability occurs in the downstream of the reconnection outflow, (2) the counter- 
streaming Weibel instability is also affected by the relativistic pressure effect, and (3) the 
Weibel instability significantly modifies the downstream reconnection structure. Since the 
Weibel instability is a micro process, it may play a role in various macro instabilities, such as 
magnetic reconnection, the RDKI, the Kelvin-Helmholtz instability, as well as collisionless 
shocks. 
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FIG. 1: (Color online) (a) Snapshot of the right half of the main simulation domain. Magnetic 
field lines (contour), plasma density (color contour) and plasma flow (arrows), (b) Out-of-plane 
field structure (B y ) in the selected region and the electric current system (arrows) in the xz plane, 
(c) Charge distribution p = [d p — d e ]/do in color. The dotted line shows E z = and the solid lines 
are contour of electric field E z with AE Z = O.IBq. 
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FIG. 2: (Color online) Time development of out-of-plane magnetic field B y along the neutral plane 
(z = 0). Profiles at three stages (i/r c = 68, 74, 80) are indicated by solid color lines. 
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FIG. 3: (Color online) Plasma distribution function in the Weibel active region (22 < x/X < 
26, — 2 < z/X < 2) at t/r c = 80 in the x-z four-velocity space, normalized by c. 
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FIG. 4: Dispersion relation of the two-dimensional purely-growing mode for k z = 5k x . The growth 
rate for three counter-streaming velocities (v = 0.6,0.65,0.7c), the cold-beam limit counterpart for 
v = 0.65c (dashed line), and observed rate (white square) are presented. 
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FIG. 5: (a) Growth rate of the two-dimensional purely-growing mode as a function of k = (k x , k z ). 
(b) The same, but for the cold-beam limit case of p = 0. 
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FIG. 6: (a-d) Field properties along the neutral plane (z = 0); B z /Bq {thick line), B y /Bo (thin 
line), and E v /Bq (dashed line) are presented. The dotted line shows B' z /Bq, obtained from a 
simulation without B y . (e-g) Properties of selected particles; (e) energy, (f) velocity, and (g) fields 
at its position are shown. 
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FIG. 7: (Color online) (a) Spatial distribution of selected positrons, whose x-velocity is slower 



than TD. (b) Spatial distribution of positrons, faster than TD. 
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FIG. 8: (Color online) (a-c) Same as Fig. 1, but for t/r c = 120 in run A [upper half) and in run 
B {buttom half), (d) x-velocity and (e) y-current at t/r c = 120 in run A (upper half) and in run 
B (buttom half). 
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FIG. 9: Energy spectra in the right half main simulation domain. 
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